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As a test of general relativity on cosmological scales, we measure the 7 parameter for the 
growth rate of density perturbations using the redshift-space distortion of the luminous red 
galaxies in the Sloan Digital Sky Survey (SDSS). Assuming the cosmological constant model, 
which matches the results of the WMAP experiment, we find 7 = 0.62 + 1.8(o"8 — 0.8) ±0.11 
at 1-sigma confidence level, which is consistent with the prediction of general relativity, 
O |. 7 ~ 0.55 ~ 0.56. Rather high value of as(> 0.87) is required to be consistent with the 

^ ' prediction of the cosmological DGP model, 7 ~ 0.68. 

h. 

Introduction - Modified gravity models, e.g., f(R) gravity, TeVeS theory, DGP 
model, have been proposed as possible alternatives to the dark energy model. Al- 
though these models might not be complete, they pose ambitious challenges to the 

■ fundamental physics, suggesting one to go beyond the standard model. In fact, a lot 
&\ \ of dark energy surveys are already under progress Testing the theory of gravity 

on cosmological scales will surely become one of the important objectives of these 

■ future large surveysPSM 

, Measurement of the growth of density perturbations will be the key for testing 

the gravity theory.0®®'^M£!IK[I3ni£!2} Several authors have already inves- 
tigated the growth of density perturbations as a way of constraining these theo- 
ries.^ JIBED j n the future weak lensing statistics will be a promising probe of the 
density perturbations, while the redshift-space distortions may also be useful for 
constraining the growth rate of perturbations PS'EB Recently, Guzzo, et al. have 
reported a constraint on the growth rate by evaluating the anisotropic correlation 
function of the galaxy sample from the VIMOS-VLT Deep Survey (VVDS)P3 

The characteristic redshift of the VVDS galaxy sample is rather large. However, 
the survey area of the VVDS sample is small. This is a disadvantage in detecting the 
linear redshift-space distortions. In the present paper, we use the Sloan Digital Sky 
Survey (SDSS) luminous red galaxy (LRG) sample from the Data Release 5, whose 
survey area is around 5000 square degrees. In this letter, we present the results of 
the multipole power spectrum analysis for the SDSS LRG sample, and subsequently 
use it to measure the 7 parameter for the growth rate of density perturbations. This 
gives a simple test of general relativity. Throughout this paper, we use the units 
where the light speed equals 1. 

Measurement of the quadrupole spectrum - The (linear) growth rate is defined by 
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f = d\nDi/dlna, where D\ is the growth factor, and a is the scale factor. Due to the 
continuity equation of the matter density the linear velocity field is related to the time 
derivative of the matter density contrast, which itself is proportional to the growth 
factor D\. The peculiar velocity of galaxies contaminates the observed redshift, 
which leads to the difference in the radial position if the redshift is taken as the 
indicator of the distance. This causes the difference in the spatial clustering between 
redshift space and real space, which is called the redshift-space distortion. The 
power spectrum including the redshift-space distortion can be modeled as (e.g.. |22|) ) 
P(k, fx) = (b(k) + f fj, 2 ) 2 P mass (k)V(k, fx), where n is the directional cosine between the 
line of sight direction and the wave number vector, b(k) is the bias factor, P mass (k) 
is the mass power spectrum, T>(k, (a) describes the damping factor due to the finger 
of God effect. 

Thus, the redshift-space distortion causes the anisotropy of the clustering am- 
plitude depending on \x. The multipole power spectra are defined by the coefficients 
of the multipole expansion^ JH* P(k,/j,) = E^=o,2,- Pe(k)C e (fi) (2£ + l), where C t (ji) 
are the Legendre polynomials F^l The monopole Po(k) represents the angular aver- 
aged power spectrum and is usually what we mean by the power spectrum. P2(k) 
is the quadrupole spectrum, which gives the leading anisotropic contribution. The 
usefulness of the quadrupole spectrum for the dark energy is discussed in [25]) . 

Within the linear theory of density perturbations, the ratio of the quadrupole 
to the monopole is given by P 2 (k)/P (k) = [4/3/3 + 4/3 2 /7]/5[l + 2/3/3 + (3 2 /5], 
where (3 = f/b(k). Thus, we can measure the growth rate from the quadrupole- 
monopole ratio. However, there are two difficulties to perform this method. The first 
is the clustering bias b(k), for which we need other independent information. The 
second is the effect of nonlinear velocity, the finger of God effect, which significantly 
contaminates the quadrupole spectrum, as one can see in Figure 1-b. 

In the present work we measured the monopole and quadrupole power spectra 
in the clustering of the SDSS DR5 luminous red galaxy sample. The galaxy sample 
used in our analysis consists of 64867 galaxies over the survey area of 4780 deg 2 and 
redshift range 0.16 < z < 0.47P^ We have excluded the southern survey stripes 
since these just increase the sidelobes of the survey window without adding much of 
the extra volume. We have also removed some minor parts of the LRG sample to 
obtain more continuous and smooth chunk of volume. 

The scheme to measure the monopole and the quadrupole spectrum is the same 
as the one described in reference 24). In this analysis we adopted the comoving 
distance of the fiducial model, which we assumed to be a flat universe with the cos- 
mological constant and matter density parameter f2 m = 0.3. In our power spectrum 
analysis we fixed the weight function equal to 1. 

Figure 1-a shows the monopole spectrum. The error bars correspond to the 
1-sigma errors. The result is consistent with the previous measurement by Hiitsi in 
reference I26p . which is also plotted in this figure for comparison. As discussed in 
reference I26p. the monopole spectrum reveals the baryon acoustic oscillation feature. 



*' Note that our definition of the multipole spectrum Piik) is different from the conventional 
definition by the factor 21+1. 
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Figure 1-b plots the ratio of the quadrupole to the monopole. The quadrupole 
spectrum changes the signature at the wave number k ~ 0.3 /iMpc -1 . This is because 
the finger of God effect becomes significant on small scales. The solid curve is the 
theoretical curve for the ylCDM model with Q m = 0.3, n s = 0.96 (initial spectral 
index), eg = 0.8, h = 0.7, 7 = 0.56 and a v = 360km/s (see also below), and assumes 
a clustering bias b(k) = 1.2[1 + 0.2(A;/0.1/iMpc" 1 ) 1 / 2 ]. 
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Fig. 1. (a, Left) Monopole power spectrum as a function of the wave number. The dark (black) 
points represent this work, the light (green) points the previous result (b, Right) Quadrupole 
power spectrum divided by the monopole spectrum. The solid curve is the theoretical curve of 
the /1CDM model. 



Testing the growth rate - By using the quadrupole spectrum, we perform a 
simple test of the gravity theory. In particular, we focus on the 7 parameter, which 
is introduced to parameterise the growth rate as, / = d\n D\{a) / d\n a = l? m (a) 7 , 
where Q m {a) = Hq J? m a -3 / H (a) 2 , H(a) = a/ a is the Hubble expansion rate, Hq(= 
100/ikm/s/Mpc) is the Hubble parameter. 

Measurement of 7 provides a simple test of the gravity theory. Within general 
relativity, even with the dark energy component, 7 takes the value around 7 ~ 0.55.^ 
However, 7 may take different values in modified gravity models. For example, 
7 ~ 0.68, in the cosmological DGP model including a self-acceleration mechanism. 
Thus, the measurement of 7 is a simple test of general relativity. 

As mentioned in the previous section, we need to take the clustering bias and 
the finger of God effect into account. For the finger of God effect we adopt the fol- 
lowing form of T>{k, //), the damping due to the nonlinear random velocity, T>(k, fj,) = 
1/[1 + (kfj,a v / Hq) 2 /2], where a v is the one dimensional pairwise velocity dispersion 
(e.g., 29)). This form of the damping assumes an exponential distribution func- 
tion for the pairwise peculiar velocity. In order to determine the clustering bias, we 
simply fix the value of erg. If <7g is fixed, and the cosmological parameters and the 
clustering bias are given, we can compute the monopole spectrum Pg heor (fc) with 
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Pg(k) = \ J^j dfiP(k, n)Ci(p), where we use the Peacock and Dodds formula for the 
mass power spectrum P mass (k)W^ We determine the clustering bias b(ki) through 
the condition PQ bs (fcj) = PQ hcov (ki) using a numerical method. Here -Pg bs (/cj) is the 
measured value of the monopole at wave number hi, and Pg heor (fcj) is the correspond- 
ing theoretical value. Figure 2 exemplifies the bias obtained by the above method 
for the cases of as fixed as as = 0.9, 0.8, 0.7. 

We used the monopole 
spectrum to determine the bias, 
and the quadrupole spectrum 
to obtain constraints on 7 and 
a v . Since the galaxy sam- 
ple covers rather broad redshift 
range, 0.16 < z < 0.47, the ef- 
fect of the time-evolution (light- 
cone effect) should be consid- 
ered proper lypQ) However, for 
simplicity, we here evaluated 
the theoretical spectra at the 
mean redshift of z = 0.31. 

Figure 3-a demonstrates 
the contours of Ax 2 in the 
7 versus a v parameter plane. 



3.0 



2.5 



2.0 



1.5 



1.0L 
0.0 




(78=0.9 



(78=0.8 



_» (78=0.7 

I I . 



0.1 



0.2 
k[hMpc" 



0.3 



0.4 



Fig. 2. The bias b(ki) obtained from our numerical 
method. The three curves correspond to as = 
0.9, 0.8, 0.7 and the other parameters are the same 
as those of Figure 1-b. For the case erg = 0.8 the er- 
rors are estimated from those of P(ki). The dashed 
curve is b(k) = 1.2[1 + 0.2(fc/0.1/iMp C - 1 ) 1/2 ]. 



We compute 



X 

T,i[P2 bs (k) - P^{ki)} 2 /[AP^ 



as x 



and AP$ bs (ki 



2 \ n i 
where P 2 obs (fc 
are the measured values and er- 



rors as plotted in Figure 1-b. 
P2 T (ki) are the correspond- 
ing theoretical values. The solid 
curves assume as = 0.8, while the dotted (dashed) ones a$ = 0.7 (as = 0.9). 
The other parameters are fixed as fi m = 0.28, n s = 0.96, h = 0.7.E3 In Fig- 
ure 3-a we plot the contour levels of Ax 2 = 2.3 (inner curves) and 6.2 (outer 
curves), which correspond to the 1-sigma and 2-sigma confidence levels of the x 2 
distribution. Clearly, the higher value of as favours higher value for 7. We find 
7 = 0.62 + 1.8(ff 8 -0.8) ±0.11 (±0.19) and a v = 367 + 80(a 8 -0.8) ± 16 (±27) km/s 
at 68 (90) percent confidence level, respectively. The value of 7 is consistent with 
general relativity. The result is not sensitive to the inclusion of the baryon oscillation 
in the theoretical power spectrum. 

The relation of 7 and as can be understood as the degeneracy between as and 
the growth rate / in the following way. As the observed power spectra can be roughly 
written as P obs oc b 2 (k)ajD 2 (z)/D 2 (z = 0) and P 2 obs oc b(k)falD 2 (z)/D 2 (z = 0), 



*' As the observed power spectra are obtained by adopting the distance-redshift relation of the 
fiducial model, the flat /1CDM model with f? m = 0.3, the cosmological distortion effect is taken into 
account properly in our theoretical computationpS'^ 
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the degeneracy between as and the growth rate / (or 7) in our method is given by 
fa%D\(z)/Di{z = 0) = constant. 

Figure 3-b is the analogue of Figure 3-a, with the expansion history now taken to 
be that of the spatially flat DGP model, which follows H 2 (a) — H(a)/r c = 8irGp/3, 
where p is the matter density and r c = Hq/(1 — Q m ) is the crossover scale related 
to the 5-dimensional Planck mass. The expansion history in this model can be well 
approximated by the dark energy model with the equation of state parameter w(a) = 
wo+w a (l—a), where wq = —0.78 and w a = 0.32, as long as Q m ~ 0.3.^ However, the 
Poisson equation is modified, and the growth history is approximated by the formula 
with 7 ~ 0.68. In order to be consistent with 7 = 0.68, Figure 3-b requires higher 
value of a 8 as compared to the ylCDM case. We find 7 = 0.47 + 1.7(er 8 - 0.8) ± 0.09 
at 68 percent confidence level, which requires a% > 0.87. 

The pair-wise velocity dispersion measured in this work is somewhat smaller than 
the theoretical model in the referenced Li, et al. investigated the pair-wise velocity 
of the SDSS galaxiesPS Their analysis is limited since it is based on galaxies with 
redshifts less than 0.3, however, they report the dependence of the pairwise velocity 
dispersion on galaxy properties and also on scales. 

Conclusion - In summary, we measured the monopole and quadrupole spectra 
in the spatial clustering of the SDSS LRG galaxy sample from DR5. The monopole 
spectrum is consistent with the previous result by Hiitsi. Using the quadrupole 
spectrum, we measured the 7 parameter for the linear growth rate and the pair-wise 
peculiar velocity dispersion. The measurement of 7 provides a simple test of general 
relativity. The measured value of 7 is consistent with general relativity, however, it 
is inconsistent with the cosmological DGP model, 7 ~ 0.68, as long as as < 0.87. 
If a constraint on as from other independent sources, e.g., the cosmic microwave 
background anisotropies, is included, we would be able to obtain tighter constraint 
on the DGP modelP^ The constraint on 7 can be applied to other modified gravity 
models, given that the value of 7 which characterizes a particular model is found, as 
discussed by Linder and Calm® (cfP^). In this work we only considered a spatially 
flat model and fixed the cosmological parameters so as to match with the results 
of the WMAP experiment. The constraint on 7 will be weakened by including the 
uncertainties of the cosmological model, e.g., the dark energy properties. 

Acknowledgements We would like to thank H. Nishioka, Y. Kojima and Y. Suto 
for useful discussions and comments. We thank G. Nakamura, K. Koyama and 
R. C. Nichol for useful comments and discussions on the earlier version of the 
manuscript. This work was supported by a Grant-in-Aid for Scientific research 
of Japanese Ministry of Education, Culture, Sports, Science and Technology (No. 
18540277). 

References 

1) A. Albrecht, et al, arXiv:astro-ph/0609591 

2) J. A. Peacock, et al, |arXiv:astro -ph/0610906 

3) D. Huterer, E. V. Linder, Phys. Rev. D 75 (2007), 023519 

4) M. Ishak, A. Upadhye, D. N. Spergel, Phys. Rev. D 74 (2006), 043513 

5) K.Yamamoto, B.A.Bassett, R.C.Nichol, Y.Suto, K.Yahata Phys. Rev. D 74 (2006), 063525 



6 



K. Yamamoto, T. Sato, and G. Hiitsi 



500 



450 



400 - 



350 - 



300 




450 



400 



350 



300 
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